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We  explore  and  characterize  the  behavior  of  ensembles  of  nitrogen  vacancy  qubits  under  transverse 
magnetic  fields.  In  this  regime,  we  predict  and  experimentally  observe  crossings  and  anti-crossings 
of  the  hyperfine  splitting  levels  which  occur  at  an  order  of  magnitude  less  field  than  the  axially- 
applied  510  G  and  1000  G  level  anti-crossing  characteristic  of  the  excited  and  ground  states.  This 
system  is  a  good  candidate  for  electric  field  sensing,  temperature  sensing,  and  quantum  memories 
in  this  novel  transverse  field  regime.  In  addition,  we  quantify  both  advantages  and  limitations  of 
addressing  distributed  ensembles  of  nitrogen  vacancies  in  a  bulk  diamond  under  transverse  magnetic 
field. 


Nitrogen  vacancy  defect  centers  in  diamond  (NVs),  are 
optically  polarizable  quantum  systems  with  high-contrast 
spin-dependent  fluorescence.  Using  electron  spin  reso¬ 
nance  under  ambient  conditions,  sensitivity  to  electric 
fields  [1,  2],  transverse  and  axial  magnetic  fields  [3-8], 
temperature  [9-12],  strain  [13]  and  pressure  [14]  have 
been  observed  via  resonance  frequency  shifts  of  the  NV 
ground-state  manifold.  When  used  as  a  sensor,  non- 
separable  sensitivity  to  multiple  fields  poses  limitations. 
However,  the  Hamiltonian  which  governs  the  energy 
shifts  can  be  tailored  to  enhance  (or  suppress)  sensitiv¬ 
ity  to  these  external  fields.  Often  touted  for  sensitiv¬ 
ity  to  magnetic  fields,  NV  energy  level  shifts  are  merely 
quadratic  ( i.e .  zero  to  first  order)  in  the  strain-dominant, 
low  magnetic  field  regime.  An  axial  or  transverse  mag¬ 
netic  bias  field  (B\\  or  B±)  lifts  this  degeneracy,  ener¬ 
getically  separating  the  spin  states  and  increasing  sensi¬ 
tivity  to  magnetic  or  electric  fields,  respectively.  In  this 
manuscript,  we  highlight  an  unexplored  low- field  regime 
where  energy  level  splittings  are  on  par  with  frequency 
shifts  of  the  applied  field.  Here  we  account  for  both  the 
electron  and  the  nuclear  spin  of  the  ensemble,  which  re¬ 
veals  complex  dynamics  of  nuclear  spin  state  degeneracy 
and  novel  avoided  level  crossings.  Note  these  features 
occur  at  low  field  as  compared  to  the  B\\  =  510  G  [15] 
and  B\\  =  1000G  [16]  excited  and  ground  state  cross¬ 
ings,  which  have  been  used  for  nuclear  spin  polarization 
to  decrease  resonance  line  width  and  increase  resonance 
contrast  for  increased  sensitivity  to  shifts  [17].  With 
exquisite  agreement  between  theory  and  experiment,  we 
paramaterize  the  regime  where  nuclear  spin  degeneracy 
and  crossings  occur  and  show  utility  toward  sensing  and 
quantum  information  applications. 

In  this  paper,  we  examine  a  transverse-mangetic-field- 
dominant  regime  where  B±  <  100  G  in  the  presence  of 
small  axial  fields  (<  3G).  The  observed  resonance  fre- 


FIG.  1.  Energy  levels  of  the  NV  system  with  B±  =  35  G 
and  B±  =  45  G  assuming  negligible  strain  and  electric  fields. 
The  non-dressed  electronic  spin  states  are  indicated  with  solid 
lines  (ms  =  +1),  dashed  lines  (ms  —  —1),  and  arrows  (ms  — 
0),  while  the  non-dressed  nuclear  spin  states  are  represented 
by  the  colors  red  (mi  =  0),  blue  (mi  —  —1),  and  yellow 
(mi  =  +1).  Note  mi  =  0  energy  levels  overlap  the  mi  =  d=l 
energy  levels  at  lower  transverse  field  as  in  (a),  but  not  at 
higher  fields  (b). 
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magnetic,  electric,  and  strain  fields  [18],  taking  into  ac¬ 
count  the  zero-field  splitting,  nuclear  and  electronic  zee- 
man  shifts,  stark  shifts,  hyperfine  splitting,  and  nuclear 
quadrupole  effects: 

Ugs  ={Dga  +  d{lUz)(S2z  -  ±S{S  +  1)) 

-  d±[nx(S2x  -S2y+  ILy(SxSy  +  SySx)\ 

+  p[ge(S  •  B)  -  gn(B  •  /)] 

+  A\\SZIZ  +  A±{SXIX  +  Syly)  +  PgS(I2  -  j) 

(1) 

where  Dgs  is  the  temperature-dependent  ground  state 
crystal  field  splitting,  d±  and  d\\  are  the  components  of 
the  ground  state  electric  dipole  moment,  the  total  effec¬ 
tive  electric  field  ft  =  E  ±  &  encompasses  both  static 
electric  and  strain  fields,  ge  and  gn  are  the  electric  and 
nuclear  Lande  g- factors,  (3  oc  fiB  is  the  Bohr  magneton 
constant,  B  is  the  applied  magnetic  field,  A\\  and  A± 
are  the  axial  and  non-axial  magnetic  hyperfine  tensors 
describing  interactions  with  local  nuclear  spins,  Pgs  is 
the  nuclear  electric  quadrupole  parameter,  and  Sk  is  the 
spin  projection  onto  the  k  axis  with  eigenvalues  0,  ±ft, 
corresponding  to  ms  =  0,±1  [18].  Figure  1  shows  the 
calculated  energy  levels  as  a  function  of  axial  field  for  a 
fixed  transverse  field. 

Solving  for  the  eigenvectors  of  this  hamiltonian  indi¬ 
cates  that  transverse  magnetic  fields  induce  state  mixing. 
This  state-mixing  affects  which  transitions  are  allowed. 
For  clarity,  we  neglect  intrisnic  strain  as  it  is  expected 
to  be  relatively  low  in  type  Ha  diamond  and  plot  energy 
levels  as  a  function  of  axial  field  (Fig.  1).  There  are 
several  interesting  features  to  note.  First,  off-diagonal 
terms  in  the  Hamiltonian  cause  the  ms  =  ±1  states  to 
mix  and  cross  each  other  when  we  seep  a  small  applied 
axial  field.  This  double  electron  spin  flip  anti-crossing  oc¬ 
curs  specifically  as  ms  =4-1  states  cross  ms  =  —  1  states 
with  identical  nuclear  spin  [18]. 

Additionally,  for  transverse  fields  under  ^  40  G,  as 
show  in  Fig.  1(a),  an  additional  avoided  level  crossing 
arises.  The  mj  =  0  levels  cross  the  ra/  =  ±1  levels, 
resulting  in  further  state  mixing  and  an  as-of-yet  unde¬ 
scribed  electron- nuclear  spin  flip  anti-crossing.  In  Fig. 
1,  state- mixing-induced  anti-crossings  are  seen  when  the 
black  lines  do  not  follow  the  solid  and  dashed  colored 
lines,  which  represent  the  undressesd  spin  states.  The 
strength  of  these  avoided  level  crossing  can  be  character¬ 
ized  by  the  deviation  of  the  energy  level  from  the  bare- 
basis  energy.  Figure  2  depicts  the  strength  of  these  split¬ 
tings  through  a  range  of  transverse  fields. 

As  predicted  in  Ref.  [18],  the  strength  of  the  dou¬ 
ble  electron  spin  flip  anti-crossing,  marked  “E”  in  this 
figure,  is  proportional  to  B\.  These  electron-spin  anti¬ 
crossings  are  centered  around  B\\  =  0,  ±0.75  G.  While 
this  anti-crossing  continues  to  increase  in  strength  as  B± 
increases,  the  second  anti-crossing,  indicated  “N”  in  this 
figure,  disappears  around  40  G  when  the  energy  levels 
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FIG.  2.  (a)  Anticrossing  strength  with  increasing  B±.  (b) 
Magnitude  of  small  axial  B  fields  at  which  the  centers  of  these 
anti-crossings  are  observed.  In  both  plots,  the  double  electron 
spin  flip  anti-crossing  is  marked  “E”  and  the  electron-nuclear 
spin  flip  anti-crossing  is  marked  “N”.  At  fields  greater  than 
B±  ~  40  G,  anti-crossing  “E”  continues  to  increase  quadrat- 
ically  and  anti-crossing  “N”  goes  to  zero  as  the  energy  levels 
no  longer  overlap. 


no  longer  overlap.  The  axial  field  about  which  this  anti¬ 
crossing  is  centered  decreases  from  roughly  ±0.5  G  at 
B±  =  15  G  to  0  G  around  B±  ~  40  G.  This  is  a  direct 
result  of  the  position  of  the  raj  =  0  energy  levels  moving 
across  the  raj  =  ±1  energy  levels.  The  strength  of  this 
electron-nuclear  spin  flip  anti-crossing  peaks  around  30 
G. 

While  the  Hamiltonian  describes  a  single  nitrogen  va¬ 
cancy,  here  we  highlight  the  experimental-theory  agree¬ 
ment  for  an  ensemble  of  up  to  101 2  nitrogen  vacancy  cen¬ 
ters  equally  distributed  along  the  four  diamond  cyrstalo- 
graphic  axes  using  a  light-trapping  diamond  waveguide 
geometry  [19].  We  lift  the  ms  =  ±1  degeneracy  by  ap¬ 
plying  a  static  magnetic  field  and  tune  the  microwave  ex¬ 
citation  to  resonantly  drive  the  ms  =  0  to  ms  =  ±1  spin 
transitions,  causing  a  drop  in  the  observed  fluorescence. 
We  assume  that  the  intrinsic  strain  across  the  sample 
does  not  change  over  the  course  of  the  measurement  and 
that  observed  shifts  in  the  zero  field  splitting  (ZFS)  are 
caused  by  temperature  or  electric  field  changes.  To  apply 
transverse  and  axial  magnetic  fields,  a  60  G  permanent 
magnetic  field  (magnitude  uniformity  of  0.75%  over  a  2 
cm3  volume  at  the  center)  from  a  Halbach  array  is  posi¬ 
tioned  around  the  sample.  In  addition,  up  to  ±30  G  can 
be  applied  in  ±X,  ±Y,  or  ±Z  (lab  frame)  using  Helmholtz 
coils.  The  four  possible  orientations  of  NV  centers  in  the 
diamond  lattice  result  in  four  sub-ensembles  experienc¬ 
ing  different  combinations  of  transverse  and  axial  fields. 
The  high  signal  to  noise  ratio  of  this  sample  allows  us 
to  observe  all  four  orientations  in  order  to  confirm  the 
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FIG.  3.  Experimental  setup.  The  diamond  sample  is  lo¬ 
cated  at  the  center  of  the  concentric  magnetic  field  sources 
and  is  provided  with  free-space  green  laser  excitation  and  RF 
excitation.  The  fluorescence  is  collected  with  a  silicon  pho¬ 
todetector.  The  four  orientations  that  make  up  the  possible 
sub-ensembles  of  NV  centers  are  shown  with  an  applied  mag¬ 
netic  field.  An  avoided  crossing  of  ~  2.6  MHz  is  seen  in 
sub-ensemble  1  when  the  perpendicular  projection  of  the  B 
field  (B i,_l)  is  equal  to  30  G  and  the  axial  projection  of  the 
B  field  onto  the  axis  (Bit\\)  is  varied  around  zero  field. 


B\\  (G) 

FIG.  4.  Plots  of  transition  energies  as  a  function  of  axial  field 
under  fixed  transverse  field.  Transitions  for  B±  =  35  G  and 
B±  =  45  G  are  shown,  which  correspond  to  the  energy  levels 
in  Fig.  1.  Hyperfine  interactions  with  nearby  13C  atoms  are 
also  visible. 


orientation  of  the  applied  magnetic  fields.  We  begin  by 
aligning  our  total  magnetic  field  to  be  perpendicular  to 
one  of  the  four  orientations  (Fig.  3)  using  the  Helmholtz 
coils.  Next,  we  apply  a  small  axial  field  sweep  in  order  to 
observe  the  anti-crossings  in  the  electron  spin  resonance 
plots. 

We  examine  the  differences  between  the  energy  levels, 
or  transition  energies,  which  we  can  view  using  optically 
detected  magnetic  resonance  (ODMR).  Figure  4  overlays 
theory  found  by  applying  selection  rules  to  the  transitions 
seen  in  the  energy  levels  in  Fig.  1  over  experimental  data. 
The  diameter  of  the  circle  corresponds  to  the  probability 
of  the  transition.  Two  values  of  applied  transverse  field 
are  shown,  corresponding  to  the  two  values  in  Fig.  1.  For 
both  B±  =  35  G  and  B±  =  45  G,  we  see  clear  double¬ 
electron  spin  flip  anti-crossings  centered  at  B\\  =  0,  ±0.75 
G  (Fig.  2b)  with  anti-crossing  strengths  corresponding 
to  those  expected  from  Fig.  2a.  As  expected,  we  also  see 


the  zero-field  splitting  resonance  increase  with  increased 
applied  transverse  field. 

Concentrating  on  the  region  of  interest  at  30  G  where 
the  electron- nuclear  spin  flip  anti-crossing  is  maximized, 
we  increase  the  resolution  of  our  scan  and  decrease  the 
RF  power  to  reduce  power-broadening  in  our  line  widths 
(Fig.  5).  The  decreased  RF  power  also  results  in  a  de¬ 
crease  in  contrast  and  consequently  a  lower  signal-to- 
noise  ratio.  However,  we  clearly  observe  the  electron- 
nuclear  spin  flip  anti-crossing  at  just  under  ±0.4  G  ap¬ 
plied  axial  field  as  predicted  by  theory. 

There  are  several  ways  that  this  moderate  trans¬ 
verse  magnetic  field  regime  and  the  anti-crossings  seen 
therein  can  impact  diamond-based  sensing.  First,  the 
double-electron  spin  anti-crossing  suppresses  sensitivity 
to  changes  in  axial  magnetic  fields,  which  addresses  an 
ongoing  challenge  in  diamond-based  sensing:  decoupling 
the  NV’s  sensitivity  to  its  many  environmental  influences. 
As  the  applied  transverse  field  is  increased,  the  magni- 
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FIG.  5.  Detail  of  double-electron  and  electron- nuclear  spin 
anti-crossings  at  B±  =  30  G.  The  RF  excitation  power  was 
reduced  by  10  dB  to  reduce  power  broadening.  This  allows  the 
higher  resolution  features  to  be  observed,  however  it  is  at  the 
expense  of  signal-to- noise-ratio.  Double-electron  spin  anti¬ 
crossings  of  magnitude  ~  2.6  MHz  are  seen  centered  around 
B\\  —  0,±0.75  G  and  electron- nuclear  spin  anti-crossings  of 
magnitude  ~  250  kHz  are  seen  centered  around  B\\  —  ±0.35 
G. 


tude  of  this  anti-crossing  increases  quadratically,  increas¬ 
ing  the  magnitude  of  axial  field  to  which  the  NV  ensemble 
is  insensitive. 

Next,  the  double-electron  spin  anti-crossing  causes  the 
hyperfine  transition  levels  to  cross;  at  these  crossings, 
we  see  increased  signal  contrast.  In  CW  diamond-based 
sensing  applications,  the  resonance  contrast  affects  sen¬ 
sitivity  linearly:  doubling  contrast  doubles  sensitivity. 
Furthermore,  as  this  anti-crossing  increases  and  the  hy¬ 
perfine  transition  levels  are  pressed  closer  together,  all 
three  levels  can  be  effectively  degenerate,  providing  a  fac¬ 
tor  of  three  increase  in  sensitivity. 

Additionally,  in  numerical  simulations,  both  anti¬ 
crossings  show  a  dependence  on  strain  and  electric  field, 


paving  the  way  for  CW  diamond-based  strain  and  elec¬ 
tric  field  sensing  based  on  tracking  the  magnitude  of  this 
anti-crossing.  This  scheme  would  rely  on  being  able  to  re¬ 
solve  the  hyperfine  transitions,  as  well  as  the  high  signal- 
to-noise  ratio  of  the  light-trapping  diamond  waveguide 
geometry.  Furthermore,  a  simple  method  could  be  devel¬ 
oped,  based  on  the  application  of  a  transverse  magnetic 
field  in  a  low  electric  field  environment,  to  characterize 
the  intrinsic  strain  across  a  diamond  sample. 

Finally,  we  expect  d^s ,  the  temperature  dependence 
of  Dgs ,  to  decrease  with  increased  transverse  magnetic 
field  [12].  Combining  this  advantage  with  the  aforemen¬ 
tioned  increase  in  contrast  results  in  a  regime  that  is  well- 
suited  for  temperature-stabilized  measurements,  relying 
on  the  ability  to  probe  multiple  sub-ensembles  with  dif¬ 
ferent  crystallographic  projections  of  the  applied  trans¬ 
verse  magnetic  field.  Comparing  the  sub-ensembles  with 
known,  varying  dependences  on  temperature  allows  the 
temperature  dependence  term  to  fall  out,  further  decou¬ 
pling  the  diamond  from  its  environment.  This  would  be 
beneficial  for  the  sensing  of  any  other  quantity. 

The  transverse  magnetic  fields  may  help  overcome  the 
influence  of  inhomogeneous  strain  present  from  diamond 
growth  processes,  a  challenge  in  working  with  very  large 
ensembles  of  NV  centers.  Combining  this  with  insensitiv¬ 
ity  to  small  axial  magnetic  fields  and  potential  insensitiv¬ 
ity  to  temperature,  the  resultant  system  also  presents  it¬ 
self  as  a  strong  candidate  for  a  quantum  memory.  Trans¬ 
verse  and  axial  magnetic  fields  could  be  used  to  read  and 
write  to  the  system,  which  is  otherwise  isolated  from  it’s 
environment. 

In  conclusion,  we  have  introduced,  described,  and 
observed  double-electron  spin  and  electron-nuclear  spin 
anti-crossings  in  the  diamond  nitrogen  vacancy  system 
under  transverse  magnetic  fields.  These  anti-crossings 
show  potential  for  a  variety  of  sensing  applications.  Sim¬ 
ilar  increases  in  contrast  can  be  achieved  to  the  much 
higher,  axially- applied  510  G  and  1000  G  level  anti¬ 
crossings  with  much  lower  applied  transverse  fields.  Fur¬ 
thermore,  this  unique  regime  allows  selective  decoupling 
of  the  NV  from  its  environment,  increasing  it’s  appeal  as 
both  a  sensor  and  a  quantum  memory. 
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